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abstract 



We discuss how the expansion and the chiral random matrix theory (xRMT) can be 

used in the study of large- A'^c gauge theories. We first clarify the parameter region in which 
each of these two approaches is valid: while the fermion mass m is fixed in the standard large- 
Nc arguments ('t Hooft large- A'^c limit), m must be scaled appropriately with a certain negative 
power of Nc in order for the gauge theories to be described by the xRMT. Then, although these 
two limits are not compatible in general, we show that the breakdown of chiral symmetry can 
be detected by combining the large- A'^c argument and the ^RMT with some cares. As a concrete 
example, we numerically study the four dimensional SU{Nc) gauge theory with Nf = 2 heavy 
adjoint fermions, introduced as the center symmetry preserver keeping the infrared physics intact, 
on a 2^ lattice. By looking at the low-lying eigenvalues of the Dirac operator for a massless probe 
fermion in the adjoint representation, we find that the chiral symmetry is indeed broken with the 
expected breaking pattern. This result reproduces a well-known fact that the chiral symmetry is 
spontaneously broken in the pure SU {Nc) gauge theory in the large- A'^c and the large- volume limit, 
and therefore supports the validity of the combined approach. We also provide the interpretation 
of the gap and unexpected A^c-scaling, both of which are observed in the Dirac spectrum. 
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1 Introduction 



From a field theoretical point of view, strongly coupled gauge theories such as quantum chromody- 
namics (QCD) are of great interest as they have a number of nontrivial phenomena in themselves. 
Because it is difficult to study those theories analytically, effective theory approaches and the 
large-A'^c limit are often considered. For example, if one considers the low-energy limit of QCD 
(and also theories with spontaneous chiral symmetry breaking (SxSB)), the form of the low energy 
effective theory is tightly constrained by the symmetries so that one obtains the chiral Lagrangian. 
If we further go to the e-regime, where the length of the box containing the system is much smaller 
than the pion Compton length, all such theories fall into one of three universality classes, which 
are exactly described by the chiral random matrix theory (xRMT). In Ref. [1] this property was 
used to demonstrate S^SB from the first principles by confirming that the spectrum of the Dirac 
operator calculated on the lattice agrees with the xRMT prediction. Another example is the 't 
Hooft large- A'^c limit p]. In particular, the large- volume independence (the so-called Eguchi- 
Kawai (EK) volume independence) [3] and the orbifold equivalence [U |5] have recently received 
large attention in the context of the lattice Monte-Carlo simulation (see e.g. [HIEIIH])- 

Given remarkable successes of these two approaches, it is natural to consider how they can be 
combined to study various QCD-like theories. Once the combined approach has been established, 
it has several interesting applications. One of them is the search for the conformal window or the 
walking technicolor model (WTM), where the approach is used to see whether chiral symmetry 
is broken or not. A candidate of the minimal WTM is the SU{2) gauge theory with two-flavors 
of fermions in the adjoint representation [9l [10]. So far, numerical results from (large- volume) 
lattice simulations indicate that this theory is inside the conformal window [llj. This situation 
may change when the number of colors increases from two. Then, the study of the large- A'c limit 
provides additional information useful to understand the phase diagram of the conformal window 
depicted in Ref. [T2] . 

In this paper, we discuss how to utilize the xRMT and the large- A^c equivalence to study the 
large- A^c gauge theories, and provide the theoretical argument for the detection of S^SB]^ It turns 
out that some cares are required because the valid parameter regions for these two techniques are 
different in general: in the xRMT limit (see Sec. [4] for the definition), the 't Hooft large- A^c 
expansion fails due to the non-trivial A^c-dependence of the expansion coefficients. Therefore, we 
cannot use the EK volume equivalence to relate the small-volume with the large-volume lattice 
theory in the ^RMT limit. However, we argue that the chiral properties of the large- volume theory 
can still be studied from the small-lattice model with finite A^c by using the analytic continuation 
provided unbroken center symmetry. 

Bearing the above motivation in mind, as a first step, we study pure SU{Nc) gauge theories on 
a 2^ lattice, aiming at taking the large- A^c limit. The goal is to find how we should extract physics 
in the large- volume limit from small lattice simulations. We call our theory on a 2^ lattice as 
the 2^-lattice model to avoid a possible confusion with the EK model for which the 't Hooft limit 
is usually assumed. Following the argument by Eguchi and Kawai in the 't Hooft limit, we let 
our lattice theory keep the center symmetry by introducing two-flavors of heavy adjoint fermions 

*In Ref. [B] and [131114) . the authors aheady used the xRMT techniques for numerical studies of chiral symmetry 
breaking of SU{Nc) gauge theories in the quenched limit and the theories with adjoint fermions, respectively. 
However, there is a subtlety (the difference of the limiting procedures explained in Sec.[4|, which makes the conclusion 
ambiguous [14]. In this paper we clarify this point in order to fully justify the method. 
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12 But notice that the theory we deal with is essentially the pure SU{Nc) gauge theory since 
they do not play any role in the low-energy dynamics. On the gauge configurations obtained, we 
calculate the low-lying spectrum of the Dirac operator in the adjoint representation and compare 
with the xRMT prediction; an agreement between them should provide the evidence for S^SB or 
equivalently non-zero chiral condensate in the 't Hooft limit. Using the EK volume equivalence, 
therefore, we conclude that chiral symmetry of pure SU{Nc) gauge theory is indeed broken. While 
the agreement of the Dirac spectrum is confirmed as expected, the A'^c-scaling of the spectrum 
turns out to be different from the one expected from the usual 't Hooft large- A'^c counting. To 
be specific, we found that the eigenvalues scale as 1/Nc rather than at a reasonably weak 

coupling. We will present a possible explanation for this phenomenon in Sec. [5} which does not 
override the occurrence of SxSB. 

This paper is organized as follows. In Sec. [2| we review the 't Hooft large-A'^c limit and its 
properties, especially the EK volume equivalence and the orbifold equivalence. In Sec. [3] we review 
the xRMT - the definition and the relationship with the e-regime of QCD-like theories with S^SB. 
In Sec. |4] we discuss the difference of the xRMT limit and the 't Hooft limit, in which we can 
use the ^RMT technique and the EK volume equivalence, respectively. Having this difference in 
mind, we explain how the xR-MT, combined with numerical simulations on a small lattice, can 
be used to study S^SB of the large- A'c gauge theories. In Sec. [5] we present the numerical results 
of the 2^-lattice model; we first confirm that the center symmetry is unbroken in the presence 
of two-fiavors of heavy adjoint fermions, and then proceed to the analysis of the Dirac spectrum 
including comparisons with the xRMT prediction and determination of the A'^c-scaling. 



2 The 't Hooft large- iV^ limit 

Let us first consider the SU{Nc) pure Yang-Mills theory, 

Sym = [ d^x TtF^,. (1) 

The 't Hooft large-A'^c limit [2] is the large-A'^c limit in which the 't Hooft coupling constant 
A = gy^Nc and the space-time V are fixed. The energy scale under consideration (e.g. distance 
between operators, the size of the Wilson loop) is also fixed. In this limit the 1 jNc expansion has 
a natural topological structure, 

oo 

{6) = Y,c,{\V)N',^\ (2) 

9=0 

where O is a properly normalized single trace operator. In the perturbation theory, contribution 
of order Nc'^^ comes from the genus-^f nonplanar diagrams, i.e. the Feynman diagrams which 
can be drawn on the two-dimensional surface with g handles (Fig. [T]). The connected correlation 
functions of more than one operators have the same structure. Therefore the l/A'^c-expansion 

^Although we added two heavy adjoint fermions, one heavy adjoint fermion is good enough to keep the center 
symmetry [T] |13| 

® More precisely, we take the couphng at some energy scale to be fixed. For example, in the case of the lattice 
regularization, we can take the bare lattice 't Hooft coupling to be the same. Then the beta function depends only 
on A at large- A^c, and hence A remains Ac-independent at any energy scale up to a 1 /Ac-correction. 
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Figure 1: Two-dimensional surface of genus g. 



is the same as the genus expansion. In the string theory, a genus-g surface corresponds to the 
string world-sheet with g closed string loops. Actually the Feynman diagrams can naturally be 
regarded as dynamical triangulations of the two-dimensional surfaces, and l/N^ can be identified 
with the string coupling constant. In Maldacena's gauge/gravity duality conjecture [15], certain 
gauge theories are explicitly related to string theories. In the large- A'^c limit, only the genus zero 
diagrams (planar diagrams) survive, or in the string terminology, the quantum string effect is 
suppressed at large-A'c. 

Next let us consider QCD with Nf fundamental fermions. The fermionic part of the action is 
given by 

Sf ermion = ^ d'^^ V/"*"^' il'^D^ + nif) (3) 

In the 't Hooft limit, in addition to the 't Hooft coupling and the space-time volume, the fermion 
mass mj must also be fixed. Because the fermions have 0{NfNc) degrees of freedom while the 
gluons have 0{N^), every time we replace the gluon loop with the fermion loop we obtain the 
factor Nf/Nc- Hence the 1 /Ac-expansion becomes 

oo 

(o) = E .''(^' ^f)^c'' ■ iNfm\ (4) 

g,b=0 

where b is the number of fermion loops. Therefore, when A^^ is fixed, diagrams with fermion loops 
are suppressed. (They are not suppressed if one considers the adjoint fermion instead.) In the 
't Hooft large- Ac limit, the parameters A, V and mf are fixed, because otherwise a nontrivial 
A'c-dependence can appear through the coefficient Cg^i,{\,V,mf). For example, in QCD, if mj 
decreases with Ac, A^c-dependent infrared divergence appears because the pion becomes lighter 
and the standard 1/A^c-counting in the 't Hooft limit can be destroyed completely. 

In the 't Hooft large- A^c limit, various nice properties hold. Below we introduce the EK 
equivalence and the orbifold equivalence. 
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2.1 The large-A/'c equivalences in the 't Hooft Hmit 
2.1.1 The Eguchi-Kawai equivalence 

Let us consider Wilson's lattice gauge theory on the periodic lattice, 

Sw = 2N^b^Yl - ^ReTrP^,(n)] , (5) 
where ^, v run from 1 to 4, 6 = 1/A = l/(c/f.^^A^c), 

n runs through the lattice, and the plaquette 

is given by 

P^^(n) = Un,^,Un+f,,uUl^^^^Ul^. (6) 

The unitary link variables Un,fi are related to the Hermitian gauge field Afj_ by C/n,/^ = 6*""^''*^"^ 
where a is the lattice spacing. The nature of the theory is characterized by the expectation values 
of the Wilson loops 

W{C) = Tr (?7„,^?7n+^,,C/n+^+,,p • • • ) , (7) 

where C is a closed path on the lattice and the right hand side is the trace of the product of the 
link variables along C. Note that the loop C can be larger than the lattice; one can just write a 
closed loop on the infinite lattice and impose the periodic boundary condition. This theory has 
the {IjNc)^ center symmetry, which multiply a phase factor to the link variables: 

C/n,M ^ e^-^^/^^C/n,^ {k^GZ). (8) 

Let us consider the 't Hooft large- A'^c limit, in which the coupling constant b and the lattice size 
L are fixed. As long as the (Zjv^)'* center symmetry is not broken spontaneously, the expectation 
values of the Wilson loops do not depend on L. This is called the Eguchi-Kawai equivalence [3]. 
One can also introduce fermions with the periodic boundary condition; for example, for both the 
fundamental and adjoint fermions, the Dirac spectrum does not depend on L when Nf and the 
mass are fixed. Therefore one can study the large-iVc theory on the infinite lattice by using a 
small lattice, say 1^ or 2^. 

In the pure glue theory, the (ZtVc)^ center symmetry is actually broken in the weak coupling 
limit 6—7-00 with fixed L [T6] . The breakdown of the center symmetry can easily be understood 
by calculating the one-loop effective action as a function of the Wilson line phases, by assuming 
four Wilson line phases to be static and diagonal. Then the off-diagonal fluctuation of the link 
variables provides an attractive interaction between the diagonal elements, so that the eigenvalues 
of the link variables favor the same value and hence the center symmetry breaks. It is nothing 
but the deconfinement transition in a very small volume. 

Introducing the adjoint matter 

The situation changes drastically when one adds the adjoint fermions, because they provide the 
repulsive force between eigenvalues. When one massless Majorana adjoint fermion is introduced, 
the theory is roughly the dimensional reduction of the four-dimensional M = 1 pure super Yang- 
Mills theory. In the continuum theory, the forces acting on the Wilson line phases cancel to all 
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order in perturbation theory. By taking into account the nonperturbative effects, both the center- 
symmetric p/7j and center-broken [18j phases can exist. (The importance of the nonperturbative 
effect was nicely demonstrated in a related context in [l9].) Whether the center symmetry breaks 
or not on a lattice is a subtle issue which depends on the detail of the lattice regularization. If we 
add more massless adjoint fermions, the center symmetry is unbroken irrespectively of the detail 
of the lattice action |17j . Furthermore, somehow surprisingly at first glance, the center symmetry 
does not break spontaneously even with very heavy adjoint fermions, whose mass is as heavy as 
the ultraviolet cutoff scale [3 [20] ; therefore one can study the pure Yang-Mills theory by using 
the EK equivalenc^ 

2.1.2 The orbifold equivalence 

Another large-A'^c equivalence, which turned out to be deeply related to the EK equivalence, 
was discovered by Lovelace He considered the S0{2Nc), USp{2Nc) and SU{Nc) Yang-Mills 
theories and found that the Wilson loops take the same expectation values in all three theories. 
This equivalence was rediscovered and generalized in the study of the string theory [5], and a 
deeper understanding in terms of the field theory was obtained [291 EO]. Today it is called the 
orbifold equivalence. 

The general statement is as follows. Let us start with the 'parent' theory, which is SO{2Nc) 
or USp{2Nc) in the case of [1]. We identify a discrete symmetry of the parent theory, which is a 
Z2 subgroup of the gauge symmetry in this example. Then we perform the 'orbifold projection' 
by keeping only the degrees of freedom which are invariant under the discrete symmetry, so that 
the 'daughter' theory {SU{Nc)) is obtained. If the projection satisfies a certain condition, the 
correlation functions of the operators in the parent theory which are invariant under the projection 
symmetry (the Z2 subgroup) and the correlation functions of the corresponding operators in the 
daughter theory take the same values, up to a calculable combinatoric factor. In particular, the 
expectation values of the chiral condensate take the same value in the SO{2Nc), USp{2Nc) and 
SU{Nc) theories with the fundamental fermions. The EK equivalence can also be understood 
as a special example of the orbifold equivalence |17j . Other valuable applications include the 
large- A'^c QCD at finite density [8] |3T1 [32l [33] , confinement in pure Yang-Mills theory [SlISS], and 
interesting properties of non-supersymmetric daughters from supersymmetric parents |36 t 137 1 [55] . 



3 The chiral random matrix theory (xRMT) 

In this section we provide a brief review of the xRMT. For more details, see e.g. |39[ I40j. 

We consider the e-regime of QCD, where the space-time volume V = is taken such that L 
is much smaller than the pion Compton wavelength and is much larger than 1/Aqcd 141] . 

-r^«L<— . (9) 

In the e-regime, the only relevant degrees of freedom are zero momentum modes of pions. Then the 
system has a universality, that is, the dynamics is determined by the symmetry breaking pattern 



^ Another way to avoid the center symmetry breakdown is the double trace deformation |21II22| . Other variants, 
quenched [16] and twisted [23] EK models, turned out to fail actually |24| |25l [26] I27j . For a further modification of 
the twisted EK model may preserve the center symmetry, see |28| . 
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and the microscopic details of the theory do not matter. Therefore QCD can be replaced by the 
XRMT, which is a random matrix model with the same symmetry breaking pattern; in particu- 
lar, the spectrum of the Dirac operator can be calculated from the ^R-MT. Note that the same 
argument holds for other QCD-like theories as long as chiral symmetry is broken spontaneously. 
The partition function of the ^RMT is given by 

Nf 

dcDjJdetP^e-^^'*--*'*, (10) 
/=i 

where $ is a x (N + v) matrix and G is a normalization parameter. N corresponds to the 
size of the system (roughly speaking the space-time volume), and v is the topological charge. 
Correspondingly to the thermodynamic limit of QCD, N is sent to infinity. In this limit, however, 
the fermion mass ruf must be scaled so that nifN, which is (roughly speaking) identified with 
nifV, is fixed. 

Note that we can define the 't Hooft large-A^ limit (not the large-A^c limit) for the xRMT, 
in which is fixed. The limit one has to take for the comparison to QCD {nifN fixed) is not 
this 't Hooft limit. This difference is crucial when we compare the large- A'^c gauge theories and 
xRMT, as we will see in Sec. m In addit ion, it is important to notice that N is identified by the 
total degrees of freedom associated with the low-energy dynamics and the individual degrees of 
freedom, such as the volume and the number of colors, does not appear in the ^RMT explicitly. 

The ensemble and the Dirac operator T) are chosen so that the Dirac operator has the same 
symmetries as the counterparts in QCD and QCD-like theories. Depending on the universality 
classes, there are three xRMTs, which are distinguished by the Dyson index /3 = 1, /? = 2, and 
/3 = 4 



13 = 2 (e.g. QCD and SU{Nc) (iVc > 3) with the fundamental fermions): 

where $ is an x (A^ -|- u) complex matrix and nif (/ = 1,2, •• • ,Nf) are the fermion 
masses. 

13 = 1 (e.g. SU{2) and USp{2Nc) with the fundamental fermions): 

P. = ( "^ll ] , (12) 



where <I> is an A^ x (A^ -|- v) real matrix. 

/3 = 4 (e.g. SU{Nc) with the adjoint fermions and S0{2Nc) with the fundamental fermions): 

mfl ^ 
nifl 



I (13) 



where $ is a 2A^ x 2(A^ -|- u) quaternion real matrix, which can be written by using four 
A^ X (A'" -|- v) real matrices (/>o, </>!, 4>2 and (p^ as 
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Because of its simplicity, the xRMT can be solved analytically. It has been applied not just 
to the test of S^SB in the lattice simulations, but also to other important problems such as the 
QCD at a finite baryon chemical potential |43l I44| and the phase structure of the Wilson Dirac 
operator [i5lH6l li7]. 

4 Large- A^c versus xRMT 

In this section, we establish the way to use the xRMT techniques in large- iVc gauge theories. For 
concreteness, we consider the SU{Nc) lattice theory with volume V and fermion mass m. Here 
the volume V is arbitrary, although we concentrate on a small fixed volume in Sec. [5| 

4.1 Difference of the limits 

In order to compare the large- A'^c gauge theory and xRMT, we must understand the difference of 
the limits, which are required for the standard l/A'^c-counting and the universality, respectively: 

• When one compares the xRMT with the gauge theory, the matrix size of the ^RMT 
is identified with the degrees of freedom in the gauge theory which are important for the 
low-energy dynamics, A'^ ~ VN^, where the constant a depends on the representation of 
the fermion in general. (As we will see, a = 1 for the massless probe adjoint fermion of 
the 2^-lattice model. Note that it is different from the usual counting in the 't Hooft limit, 
a = 2.) In order for the universality to hold, mN ~ mVN" must be fixed as we take the 
large- A^c limit. Let us call it the '^RMT limit'. This limit is compatible with the condition 
for the e regime in Eq. [9j 

• For the standard 1/Nc counting, the ordinary 't Hooft large- A'c limit, in which m and V are 
fixed, must be taken. In this limit the large- A'^c equivalences (e.g. the EK equivalence) hold. 

The xRMT limit is different from the 't Hooft limit and the standard 't Hooft l/A'c-counting 
does not hold in this limit. In order to see this, let us consider the A:-point connected correlation 
function in the xRMT (see e.g. [i2] ) 

m)'/N)conn,RMT = (^Y.(^^y) , (15) 

\ i=l ^ * ^ / RMT 

where Aj are eigenvalues of the Dirac operator. In the 't Hooft counting, {{i^ip)^ / N) connfiMT is of 
order A^^. It is true when m is of order one: l/(Ai + m)^ is of order one and the summation over 
2A^ order-one quantities is of order N . However when m scales as l/N , the smallest of Aj + m 
is also of order 1/A^, and hence the correlation function becomes of order N^~^. This divergence 
is analogous to the infrared divergence in SU{2>) QCD in the e-regime. (Note that tpip in xRMT 
corresponds to j d^x'f/J'i/' in QCD.) In large-A'c field theories, this corresponds to the divergence 
with a certain power of N^, which is different from the 't Hooft counting. This peculiar behavior 
can also be understood in terms of the large- A^c gauge theory; because the coefficients Cg^h in @i 
are functions of m and V , they can have nontrivial A'^c-dependences in the ^RMT limit, where m 
and V are scaled with Nc- 
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C^greement in the x RMT lim 
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Large-Nc gauge theory 



Nonzero chiral condensate 
at m^O after V^infty and/or Nc 



Figure 2: The agreement between the ^RMT and the large- A'^c gauge theory in the ^RMT limit 
assures SxSB in the 't Hooft limit and/or the large- volume limit. 



Because of this difference, the large- A'^c equivalences do not hold in the xRMT limit. As an 
example, let us consider the chiral condensate in SU{Nc), SO{2Nc) and USp{2Nc) theories. In 
the planar limit, because of the orbifold equivalence (Sec. 2.1.2[ ), they are the same as a function 
of m [SI El]. On the other hand, in the xRMT limit they can be calculated by using the xRMTs 
as a function of mN, and they do behave differently [48^ I49j . It means the orbifold equivalence 
does not hold there. (Note that it is the case even in the quenched theory, which we numerically 
study in this paper.) 



4.2 The comparison 

First let us recall how one can realize SxSB in numerical simulations of QCD in the e-regime. 
The criterion for SxSB is whether the Dirac spectrum agrees with the xRMT prediction. Once 
they are found to agree, the value of the chiral condensate is determined. Notice that the value 
of the chiral condensate obtained in the e-regime is the same as the one obtained by taking the 
chiral limit after the large- volume (or thermodynamic) limit. It means that, whatever the order 
of the limiting procedures is, we can determine the chiral condensate. 

The same logic should hold for large A'^c gauge theories. But this time, the role of the volume 
in the xR-MT limit can be played by Nc- Namely, the xRMT limit is achieved by taking the large 
Nc limit with mVNc fixed. Then, the criterion for SxSB is again the agreement of the Dirac 
spectrum in the xRMT (Fig. [2]). 

In Fig. [3] we show how one can combine the xRMT and the EK equivalence to see SxSB. We 
compare the spectrum of the Dirac operator of the probe fermion with the mass m ~ A'^^" in 
the V = 2'^ lattice and that in the xRMT; the agreement between them implies the breakdown 
of chiral symmetry (nonzero chiral condensate) in the 't Hooft limit. Using the EK equivalence, 
therefore, we conclude that chiral symmetry is spontaneously broken in the large- volume lattice 
gauge theories. The constant a is determined by the effective degrees of freedom in the low mode 
region. As we will see, in the 2^-lattice model with the massless probe adjoint fermion, small 
eigenvalues scales as 1/Nc- Therefore we take a = 1 We also found a nice agreement of the 

* In [T3] the same scaling has been aheady found for the Nf = 1 theory with a very small mass (which is 
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The X RMT 



agreement in the x RMT lim 



1 



Nonzero chiral condensate 
on the 2" lattice 
at m^O after Nc 



the 2^ lattice 



Nonzero chiral condensate 
on large-volume lattice 
at m^O after Nc -^'^ 




Figure 3: The agreement between the ^RMT and the 2^-lattice model m the xRMT limit imphes 
SxSB or equivalently non-zero chiral condensate in the 't Hooft limit. Then, one can use the EK 
equivalence to conclude S^SB in the large-volume theory. 

low-lying Dirac spectrum, which we interpret as the presence of S^SB. 

Here let us discuss the Nc scaling of the chiral condensate for later use. In the large-volume 
theory, the chiral condensate E of the adjoint fermion is related with the spectral density p(X) by 
the Banks-Casher relation 1501. 



S = |(V;V>l=MO), with p(A) = ^lX^<5(A-A„)^, (16) 

where An is the Dirac eigenvalues. The spectral density at A = is proportinal to the inverse of 
spacing of the near-zero Dirac eigenvalues, A A = Aj+i — Aj. A A is expected to scale like ~ 1/A^, 
where is the number of degrees of freedom important for the low-energy dynamics. Since the 



definition of p{X) in (16) is already normalized by V, p thus defined scales as increases. In 
the 't Hooft limit, the degrees of freedom of both gauge and fermion parts are 0{N^) and as a 
consequence S = N^T, in the 't Hooft limit, where an 0{N^) quantity S is defined by the properly 
normalized operator as 

^ = ^\{trim)\=^m, with piX) = /^^5iX-Xn)\. (17) 

\ n I 

In the xRMT limit, however, we have to carefully count the number of degrees of freedom asso- 
ciated with the low-energy dynamics which may be different from that in the 't Hooft limit; in 
general the properly normalized spectral density would be related with the spacing of small Dira 

essentially massless). In that paper, however, the authors regard that the deviation from the 't Hooft counting 
a = 2 suggests the absence of SxSB. 
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The X RMT 



The / RMT 



the 24 lattice 




agreement in 
the X RIVIT limit 



the xRMT limi 



agreement in 





the large- 
volume lattice 



Nonzero chiral condensate 
on the 2" lattice 
at m^O after Nc -^'^ 




Nonzero chiral condensate 
on large-volume lattice 
at m^O after Nc 



Figure 4: Another interpretation. The agreement between the xRMT and the 2^-lattice model 
in the xRMT hmit translates into the agreement the xRMT and the large-volume lattice thanks 
to the absence of the phase transition (the center symmetry breakdown). One can conclude SxSB 
of the large-volume theory from the agreement between ^RMT and the 2^-lattice model in the 
XRMT hmit. 

eigenvalues by p(0) = l/(yA^"AA), where the a mentioned above is determined from how the 
low- lying eigenvalues scale with N^. 

Strictly speaking, there is a subtlety in the ^RMT limit of the 2^-lattice model: because the 
space-time volume of this lattice is very small, usual derivation of the xRMT from QCD via the 
chiral perturbation theory may not be applicabl^ In order to circumvent this subtlety, one can 
take another path (Fig. |4]) as follows. First let us consider a sufficiently large volume, where the 
usual derivation of the xRMT is valid when chiral symmetry is spontaneously broken. Then, 
the distribution of the Dirac eigenvalues should agree with that from the xRMT after properly 
normalizing the eigenvalues by \{V N°')Ti. Therefore, as long as V is large enough to justify the 
chiral perturbation, ^/-dependence does not appear manifestly. Now let us shrink the volume 
further. If we consider the pure Yang-Mills without the (heavy) adjoint fermion, there is a phase 
transition (the breakdown of the center symmetry), beyond which one cannot expect the same 
expression for the distribution of the Dirac eigenvalues. However in the present case, because 
there is no phase transition thanks to the heavy adjoint fermion, it is expected that the same 
expression for the eigenvalue distribution holds at very small volume as a function of \{V N^)Ti, 
where the value of a could be different from that in the large-volume theory. (In order to confirm 
this argument a quantitative study of S on a large lattice is required and therefore we leave it as 
our future work.) Then, the agreement between the xRMT and the 2^-lattice model in the RMT 
limit translates into the agreement between the xRMT and the large-volume lattice provided 

^ One might think the standard mapping rule between the planar sector of the small-volume (i.e. EK model) 
and the large- volume theory can be used. However nonplanar diagrams can contribute in the xRMT limit. 
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unbroken centery symmetry. At large volume and in the 't Hooft limit, one can conclude S^SB in 
the ordinary manner. By further using the EK equivalence, S^SB of the 2^-lattice model in the 
't Hooft limit can also be concluded. 



5 Numerical simulation of the 2^-lattice model and comparison 
to xRMT 

In this section, we apply the strategy explained in Sec. |4]to the 2^-lattice model and detect S^SB. 



After introducing the lattice action and simulation details in Sec. 5.1, in Sec. |5.2| we confirm that 



the center symmetry, which is crucial for the use of the EK equivalence, can be kept unbroken by 



adding heavy adjoint fermions. Then, in Sec. 5.3 we calculate the Dirac spectrum and compare 
with the xRMT prediction. 

For an earlier work along the same direction, see [13]. See also [51^ I52j for recent numerical 
simulations of a single-site model. 

5.1 Lattice action and simulation details 

We take the standard the Wilson plaquette gauge action Eq. [5] 

5, = 2Af2^j;^('l-^ReTYP^.(n)y n € 2\ (18) 

n /i<u V / 

In order to preserve the center symmetry, we add two- flavor of heavy adjoint fermions, for which 
we choose the plain Wilson fermion, 

^/ = E E ^n,i Un, - « E - 7M)On+M,, + (1 + l^W^^-'^tl^n-,,] ) , (19) 
j=l n \ n / 

where b and k represent the inverse of the 't Hooft coupling, b = l/(gy^jA^c)) and the hopping 
parameter, k = (2moo -|- 8)~^ (where ttiq is the bare fermion mass), respectively. The plaquette 
P^u is defined by Eq. [6j For the fermionic action, the link variables in the adjoint representation 
are defined by 

u;^f = lTV[r^c/r^c/t], (20) 

where Tp are SU{Nc) generators in the fundamental representation. This action is invariant under 
SU{Nc) local gauge transformation 

C/n,M f^nC/n,M^^n+^, ^^n G SU{N,), (21) 



as well as (Zat^) global center transformation Eq. |8j Throughout this work, we impose periodic 
boundary conditions in all directions for both link variables and fermion fields. 

Our lattice simulations consist of two parts: 1) quenched calculations (k = 0), where rriQa is 
infinite, as a nontrivial check of our numerical code by confirming the breaking of center symmetry 
at weak coupling, 2) simulations with two heavy adjoint fermions {k = 0.09), where the center 



12 



Table 1: Simulation parameters. 
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symmetry is unbroken even at weak coupling while the low-energy dynamics are essentially the 
same as those of pure SU{Nc) gauge theory since niQ is of order of the inverse lattice spacing. 

Simulation parameters are summarized in Table [TJ We performed simulations at b = 0.5 
(weak coupling) for up to Nc = 16, which is relatively smaller than that in previous single-site 
model simulations [7l|20l|51]- As we will see in Sec. 5.2.2, however, we could obtain good large- 
Nc limits since we have additional suppression of the finite volume effects thanks to the larger 
volume of a 2^ lattice. For Nc = 8 and k = 0, we also performed simulations at 6 = 0.3 and 0.4 
corresponding to the strong and intermediate couplings, respectively. In addition, we performed 
simulations at b = 0.2 (strong coupling) for k = 0.09 and Nc = 6, 8. We used the Hybrid Monte 
Carlo (HMC) algorithm for all lattice simulations with the plain leap-flog integrator, where the 
step size is tuned such that the acceptance ratio is in the range of 70 - 80 %. The simulation 
codes were developed from the one unsed in [S3] with appropriate modification to SU (Nc) adjoint 
fermions with arbitrary large Nc- After 200 trajectories for thermalization, Nconf = 138 — 600 
configurations are accumulated for each ensemble, where every two adjacent configurations are 
separated by 10 trajectories. Statistical errors are calculated by using the standard bootstrapping 
technique. 



5.2 {'^NcY Center Symmetry 

As explained in Sec. |2.1.1[ the non-trivial condition for the large- A'^c EK equivalence is that the 
(ZtVc)^ center symmetry of the reduced volume theory must be unbroken. The presence of the 
center symmetry is established by checking the followings: (1) the Polyakov loop radially scatters 
in the vicinity of origin in the complex plane, (2) the magnitude of the Polyakov loop approaches 
zero as Nc increases, (3) the average plaquette value measured in the reduced model agrees with 
that in the large-volume lattice gauge theory. In this section, we present our findings for pure 
SU (Nc) gauge theory and the theory with two heavy adjoint fermions. 
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Figure 5: Scatter plots of the Polyakov loops for quenched simulatfons with b = 0.5. The number 
of colors are Nc = 4, 8, and 16 from the left to right. 
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Figure 6: Scatter plots of the Polyakov loops for quenched simulations with Nc = 8. The values 
of the couplings are b = 0.3 and 0.4 from left to right. 



5.2.1 Pure SU{Nc) gauge theory 

of the Polyakov loops defined by 

P^=i{y,z,w) = -^TrU^=i-o^y^z,wU^=i-i,y,z,w, (22) 

and similarly for /i = 2,3 and 4, where b = 0.5 and Nc = 4,8, 16 from left to right. In this weak 
coupling regime, the plots clearly show the center symmetry breaking as the Polyakov loops are 
localized at the elements of the center of SU{Nc), 2in^iT /Nc^ where = 0, 1, • • • ,Nc — 1. For a 
given number of configurations, the number of clusters decreases to one as the number of colors Nc 
increases from 5 to 16; this means the tunneling transitions between different center-symmetry- 
broken vacua will eventually disappear at A^c — >• co- 

For Nc = 8, we performed two more simulations with smaller values of b. The results of the 
Polyakov loops are shown in Fig. [6j At 6 = 0.3, the Polyakov loops develop a cluster around the 

^"For all scatter plots, We used about a hundred ensembles and chose one of the direction /i out of four directions. 
We found similar scatter plots for three other directions. 



In Fig. pi we show scatter plots 
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Figure 7: Average plaquette values for quenched simulations with b = 0.5. The red solid line 
represents the uncorrelated fit of the data to the function, cq + ci/N^, and we obtained cq = 
0.72733(12) and ci = 0.3325(46), where the chi-square/d.o.f is 0.67. 



1.0 



0.5 



5:^ 0.0 
S 



-0.5 



-1.0, 





• 4^ 


1 • , 





















l.Or 

0.5 
Si 0.0- 

e 

-0.5 




1.0 -0.5 0.0 0.5 1.0 



-1.0, 



1.0 



0.5 



5^ 0.0 
E 



-0.5 



1.0 -0.5 0.0 0.5 1.0 
Re(P^) 



-1.0 



1.0 -0.5 0.0 0.5 1.0 
Re(P^) 



Figure 8: Scatter plots of the Polyakov loops at 6 = 0.5 with two heavy adjoint fermions. The 
number of colors are Nc = 4, 8, and 16 from left to right. 



origin, while at 6 = 0.4 they spread out and are localized at the center-symmetry-broken vacua like 
at 6 = 0.5. Therefore, we conclude that the center symmetry is restored in the strong coupling 
regime. The boundary between the strong and weak coupling regimes is located somewhere 
between b = 0.3 and b = 0.4, which is consistent with the results in |54j . 

In Fig. [7j we plotted average plaquette values for Nc up to 16. The measured plaquette values 
turn out to scale as l/N^ and thus we perform a fit to the data using a constant plus quadratic 
function of N^^^ (red solid line in the figure). We obtained 0.72733(12) in the large Nc limit, 
which is larger than 0.7182, the value from large- volume lattice gauge theory Therefore we 
reproduced the well-known fact that the large- A'^c volume reduction for pure Yang-Mills theory 
fails at weak coupling due to absence of the center symmetry. 
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Figure 9: (Left) Average plaquette values and (Right) average values of the magnitudes of the 
Polyakov loops along with fit results in Table [2] 



Table 2: Fit results for the Polyakov loop and plaquette values. 





Data set 


Fit function 


Co 


Cl 


C2 


xVd.o.f 


Plaquette 


Nc = [4, 16] 
Nc = [2, 16] 


Co + Ci/Nc + C2/N^ 

CQ + ci/Nl + C2/Nt 


0.72053(71) 
0.72067(14) 


0.003(13) 
0.405(12) 


0.383(50) 
-0.283(61) 


0.47 
0.49 


Polyakov loop 


Nc = [5, 16] 


Co + ci/Nc 


0.0026(12) 


0.819(15) 




0.89 



5.2.2 SU{Nc) gauge theory with two heavy adjoint fermions 

Now let us add two heavy adjoint fermions in order to keep the center symmetry unbroken even 
at weak coupling. At 6 = 0.5, we performed a set of simulations for various values of Nc with 
fixed value of k = 0.09 as shown in Table [TJ In Fig.[8j we show scatter plots of the Polyakov loops 
for Nc = 4, 8, and 16. The clustering of the Polyakov loops around the origin for A'^c = 8 and 16 
clearly shows that the center symmetry is intact, which is in contrast to the case without adjoint 
fermions. 

In Fig.[9| we show the average plaquette values versus N'"^ and the magnitudes of the Polyakov 
loops versus Nc^- For the plaquette values, we examined two different fits, one fitting the data 
of 4 < A'c < 16 to Co + ci/Nc + C2/A^c (blue sohd fine) and the other fitting those of 2 < A^c < 16 
to Co + c\/Nc + C2/Nc (red solid line). For the Polyakov loops, we fit the data of 5 < A'c < 16 to 
Co + ci/Nc to obtain the large A'c limit. The fit results are summarized in Table [2] 

As discussed in [20] and [51] in great details, the leading correction to plaquette in the large- 
Nc limit for a single-site lattice model is 0(1/A^c)) instead of 0{1/N^) (in the usual 't Hooft 
counting), due to the contributions of diagonal zero modes. If we apply this argument to our non- 
single site reduced model, the 1/Ac correction term is suppressed by l/V. Indeed, we obtained a 
consistent result where the one-sixteenth of the coefficient of 1 /Nc calculated in [20] is within 
the uncertainty of our results shown in Table [2] (first row). The extracted plaquette value also 
agrees with that obtained in a single-site model [20], but it is systematically larger than that 
from the large-volume lattice calculation for pure Yang-Mills. This tiny difference comes from the 
presence of heavy fermions. The magnitude of the Polyakov loop goes to zero as Nc increases; 

^'^The value of the coefficient is 0.22 ± O.Of , which was not presented in the original paper. 
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it scales as 1/Nc in the asymptotic region. Therefore, we conclude that the center symmetry is 
unbroken for a given lattice parameters and thus the EK volume equivalence is applicable. 



5.3 Comparison to ^RMT 

As discussed in Sec. |4| to detect S^SB we compare the low-lying Dirac eigenmodes of the 2^-lattice 
model with those of xRMT in the limit of A'c — )• oo with mVN^ fixed. The simplest way to 
achieve the ^RMT limit without losing the generality might be taking m = and — )• ocp^ For 
this purpose, we calculate the low- lying spectrum of the overlap-Dirac operator V for a massless 
fermion in the adjoint representation. The operator T) is defined by 



V = M 



1 + 75 ^-^-^^ 



^H^{-M)H^{-M)^ 



(23) 



where Hyj{—M) is the hermitian Wilson-Dirac operator and the parameter M is taken to be 1.6 
in the most cases. The overlap-Dirac eigenvalues lie on a circle in the complex plane |551 156j . 
To compare the Dirac spectrum with the xRMT, we consider the projection of the eigenvalues to 
the imaginary axis, 

A. = ^"^[^^J ■ (24) 
1 - Re[Afc]/mo 

Note that the Dirac eigenvalues for adjoint fermions are appearing as conjugate pairs with 2-fold 
degeneracy and we take the distinct eigenvalues on the upper-half plane for our numerical results. 
Unless otherwise noted, the numerical results in this section are restricted to the case of 6 = 0.5 
(weak coupling). 

5.3.1 Chiral symmetry breaking 



As seen in Sec. 5.2, the adjoint fermion plays a role of the center symmetry preserver in the 2^- 
lattice model. Since the EK volume equivalence is valid for the same lattice parameters such as 
the bare coupling and the fermion mass, the equivalent large-volume theory also has the fermion 
mass of order 0(l/a) and approximates the quenched theory. As a comparison, therefore, we 
consider the xRMT for the quenched theory. We restrict the xRMT predictions to the case of 
zero-topological charge. Accordingly, the configurations yielding exact zero mode(s) are omitted 
from the analysis. 

The adjoint QCD with any number of flavors belongs to the universal class of the Chiral Gaus- 
sian Sympletic Ensemble (ChGSE). However, we also consider two other universal classes, Chiral 
Gaussian Orthogonal Ensemble (ChGOE) and Chiral Gaussian Unitary Ensemble (ChGUE), in 
order to make the comparison manifest. The distributions of the lowest eigenvalue are 



P{z) 



2±.g-(./2)-(.V8) for ChGOE 

,2 ; 



e"^ for ChGUE , (25) 

fz(3/2)/3/2(z)e-^'/2 for ChGSE 



Although this Umit looks like the 't Hooft limit (m is fixed), actually one should not regard it so because the 
infrared (IR) regulator is assumed in the usual 't Hooft counting; for the usual 't Hooft counting one must introduce 
nonzero fixed m as an IR regulator and takes the large- A^'c limit, and then sends the mass to zero. 
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Figure 10: Low- lying Dirac spectrum for Nc = 8 (left) and Nc = 16 (right). The Dirac eigenvalues 
are rescaled by matching the ensemble average of the lowest eigenvalue with the expectation 
value of the lowest eigenvalue in the xRMT for the ChGSE. The solid and dotted lines represent 
the xRMT predictions of the distributions of the lowest eigenvalue and the spectral density, 
respectively. The colored histograms are for individual eigenvalues, while the black histograms 
are for all eigenvalues. 



and the spectral densities [l9] are 



f [J^{z) + J2(z)] + i Jo(z) [1 - Jq dt Jo(t)] for ChGOE 
I [J^{z) + Jf{z)] for ChGUE 

z [Ji{2z) + Jf{2z)] - ^ Jo(2z) f^' dt Jo(t) for ChGSE 



(26) 



Histograms of the lowest twelve Dirac eigenvalues for Nc = 8, 16 at 5 = 0.5 and k = 0.09 are 



shown in Fig. [10} For the lowest eigenvalue, we found that its distribution is well described by the 
xRMT for the ChGSE (solid blue curve) after introducing a rescaled eigenvalue z hj z = XVT, to 
fit the data, where V = 2'^ and E is a free parameter. As a comparison, we show the distribution 
of the lowest eigenvalue predicted by the xRMT for the ChGUE (solid red curve) and ChGOE 
(solid black curve) using the same parameter z for the ChGSE. In addition, we plot the spectral 
densities for the ChGSE (dashed blue curve), which are in good agreement with a few lowest 
eigenvalues. According to Fig. [Sj this result would be a strong evidence that chiral symmetry of 
quenched large Nc gauge theory is spontaneously broken. 

In Fig. 11, we compare 5{Xk)/{Xi) for Nc = 16 (our largest value of Nc) with the xRMT pre- 
diction, where the braket (• • • ) represents the ensemble average. The low-lying Dirac eigenvalues 
perfectly agree with the xRMT prediction, which adds a further evidence for S^SB. We have also 
studied the strong coupling region, b = 0.2, and similarly found a good agreement with the ^RMT 



prediction (see Fig. 15). However, note that this point is actually stronger coupling than the bulk 
phase transition; this phase is not related to the continuum limit a priori. 



5.3.2 A^c-scaling and gap of the Dirac eigenvalues 

In Sec. |4| we argued that the xRMT limit of the 2'^-lattice model is analogous to that of ordinary 
QCD by replacing mV with mVN^, where the exponent a can be determined so that N^ is the 
number of degrees of freedom important at low energy, which is proportional to the inverse of the 
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Figure 11: Spacing between the adjacent Dirac eigenvalues normalized by the lowest eigenvalue 
(Ai) for Nc = 16, where 6{Xk) = {Xk) — (Afc-i) and (Aq) = 0. The red dashed lines represent the 
xRMT prediction for the /3 = 4 universal class. 
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Figure 12: (Left) Spacing between the adjacent Dirac eigenvalues multiplied by Nc for Nc = 8, 12 
and 16. (Right) Spacing between the first and second Dirac eigenvalues multiplied by Nc for 



Nr 



4,6,8,12 and 16. 



eigenvalue density around the origin. Then, because the low-lying eigenvalues scale as 1/Nc as we 
will see below, we obtain a = 1. 



In the left panel of Fig. 12 , we plot the spacing between the adjacent low-lying Dirac eigenvalues 
multiplied by Nc for Nc = 8, 12 and 16. We see a nice agreement of the data points at up to A; = 2 
for Nc = 8 and at up to = 4 for Nc = 12, where the spacings are expected to show a plateau 
for k > 2. This result implies that the near-zero eigenvalues, which are expected to reproduce the 
xRMT prediction well, scale as l/Nc- (The same scaling had been reported in [Ijj, by using the 
single-site model with a very light dynamical overlap adjoint fermion.) The right figure of Fig. 12 
shows 6{X2) multiplied by Nc for various Nc, for given statistics, they agree to each other for 
Nc > 8. We find that the eigenvalue spacing deviates from the ^RMT prediction as we increase 
k or decrease Nc, the distribution has a long-tail in the direction of the large eigenvalue. This 
behavior can be understood as follows. As we will see below, the gap appears between the {Nc — 1)- 
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Figure 13: Low-lying overlap-Dirac eigenvalues in units of radian for Nc = 6, 8, 12. 



th and A'c-th eigenvalues, where the spectral density is zero. The repulsion between eigenvalues 
becomes weaker as we approach the gap, and thus the distribution develops a long-tail to the gap. 
Therefore, the eigenvalue spacing near the gap becomes larger than expected. Note also that the 
1/Nc correction takes a rather complicated form due to the peculiar distribution; for 6X2, though 
at Nc = 4, 6, 8 the corerction looks l/Nc, at > 8 this behavior disappears and the value stays 
almost constant. For JA3 we observe a similar l/A'c-like behavior at Nc = 8,12 and 16, but we 
expect the value is saturated at Nc ^ 16. 

To describe this gap clearly, we plot the lowest twelve overlap-Dirac eigenvalues, which lie on 
a circle in the complex plain, in units of radian for A'^c = 6, 8 and 12 in Fig. [T3| The eigenvalue 
abruptly jumps at A: = 5, 7 and 11 for Nc = 6, 8 and 12, respectively. (A similar gap had also 
been found in the Nf = 1 theory [2].) For k < Nc, the eigenvalue spacing roughly scales as 
1/Nc or equivalently the density scales as Nc- The gap persists in the large- A^c limit: the A'^c-th 
eigenvalues are reasonably well fit by the function 0.3 + A.GNc^- The schematic diagram of the 
overlap-Dirac eigenvalues on a circle in the complex plain is shown in Fig. 14 The A'c-dependence 
of the eigenvalue density is expected to change from Nc to N^ as k changes from k < Nc to k > Nc- 

A possible explanation of the l/A'c-scaling of Dirac eigenvalues and the appearance of a gap 
relies on the perturbative analysis around the diagonal Wilson lines for a weakly coupled gauge 
fields. In the perturbation theory, the spectrum of the theory is governed by the background of 
the Wilson lines and thus the low-lying Dirac spectrum should be closely related to the zero modes 
|20j . For the fermions in the adjoint representation the number of zero modes of the Wilson lines 
is {Nc — 1), which is the effective number of degrees of freedom, while the number of total degrees 
of freedom is N^. According to our argument in Sec.[4| therefore, the low-lying Dirac eigenvalues 
should scale as 1/Nc, being consistent with our numerical results. We emphasize that this analysis 
is only true for weakly coupled large A'' gauge theory in a compact space. 

When Nf = 1, the position of this gap has been studied for several values of the coupling 
constant 114J, which turned out to be almost independent of the coupling constant in the lattice 
unit. Therefore, in the physical unit, the scale corresponding to the location of the gap diverges 
as the lattice cutoff increases. This result looks natural, because a new physical scale appears 
otherwise. 

It is also satisfactory from the universality point of view: SU{Nc) theory with Nf adjoint 
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Figure 14: (Right) Schematic diagram of the eigenvalue density of the overlap-Dirac operator on 
a circle in the complex plane. 



fermions and SO{2Nc) theory with Nf fundamental fermions are equivalent in the ^RMT limit, 
because they are described by the same xRMT. However they are completely different in the 't 
Hooft large- A'c limit; whereas the former has 0{N^) fermion degrees of freedom, the latter has only 
0{Nc). In order for them to become identical in the xRMT limit, the fermionic degrees of freedom 
must match somehow. But now we know the mechanism: only 0{Nc) degrees of freedom survives 
in the low energy limit of the adjoint theory, because — Nc eigenvalues become infinitely large. 
(In the case of the fundamental fermions in the probe limit, there is no gap [6j). 

At strong coupling phase, the perturbative treatment around the diagonal background Wilson 
lines is no longer reliable: the zero modes can be lifted by gauge fluctuation and the off-diagonal 
components would be the same order of magnitude of the diagonal components. In contrast to the 
case of weak coupling, therefore, we expect that the gap is absent and the eigenvalue spacing is of 



order 0{1/N^). In the left panel of Fig. 15, we plot the eigenvalue spacing multiplied by for 



Nc = 6 (red circle) and Nc = 8 (blue square) and at 6 = 0.2. The data shows a nice plateau and 
agrees to each other, implying that the eigenvalue scales as l/N^ without any gap in the strong 
coupling regime, as expected. In the same manner, it is expected that the Dirac eigenvalues scale 
as 1/N^ and the gap does not exist even at weak coupling if the volume is sufficiently large: as the 
volume increases the momentum gap between the lowest and the first excited state, which differ 
by (27r/L), decreases and at some point the zero modes are lifted by gauge fluctuation. 

6 Conclusion and discussion 

In this paper we considered how to apply the xRMT techniques to large- A^c gauge theories. After 
giving general considerations, we provided a numerical demonstration by using the 2^-lattice model 
as an example. The most important lesson is that the 't Hooft large- A'c limit and the xRMT limit 
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Figure 15: (Left) Spacing of the adjacent Dirac eigenvalues multiplied by for A'c = 6, 8 with 
b = 0.2, where 6{Xk) = {Xk) — {^k~i) and (Aq) = 0. (Right) Spacing between the adjacent Dirac 
eigenvalues normalized by the lowest eigenvalue (Ai). The red dashed lines represent the ^RMT 
prediction for the /3 = 4 universal class. 



(the microscopic limit) are not compatible in general: the former is the large- A'^c limit with the 
fermion mass m and the space-time volume V fixed, while the latter requires mVN^ to be fixed, 
where a is a positive constant which may depend on the theory. The value of a is a unity in the 
example we studied in Sec. 5, which is different from the usual 't Hooft counting. 

An important consequence of the difference between two limits is that several properties in 
the 't Hooft large- A^'c limit (e.g. the equivalence between SU{Nc), SO{2Nc) and USp{2Nc) gauge 
theories) do not hold in the xRMT limit. This fact must be appreciated when one applies the 
large- A^c and/or xRMT approaches to QCD and related theories; although these two approaches 
provide valuable 'exact' results, they are valid in different parameter regions and hence one has 
to carefully choose more suitable method depending on the physics he/she studies. In spite of 
the difference of two limits, one can still detect S^SB of large-A'^c gauge theories as we have 
demonstrated in Sec. [H 

Rather curiously, we observed a nice agreement between the ^RMT and 2^-lattice model even 
when the center symmetry in the latter is broken spontaneously. Of course we cannot relate this 
fact to SxSB in the large-volume theory, because neither the EK equivalence nor the analytic 
continuation can be used. This is presumably because the space-time dimension is not important 
for the universality argument; only the pattern of SxSB matters. 

As a next step, we are studying whether the SU{Nc) gauge theory with dynamical adjoint 
fermions goes through SxSB or not, with an application to technicolor models in mind. We hope 
to report the result in near future. 
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